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We show that odd-parity superconductivity occurs in multilayer Rashba systems without requiring 
spin-triplet Cooper pairs. A pairing interaction in the spin-singlet channel stabilizes the odd-parity 
pair-density-wave (PDW) state in the magnetic field parallel to the two-dimensional conducting 
plane. It is shown that the layer-dependent Rashba spin-orbit coupling and the orbital effect play 
essential roles for the PDW state in binary and tricolor heterostructures. We demonstrate that the 
odd-parity PDW state is a symmetry-protected topological superconducting state characterized 
by the one-dimensional winding number in the symmetry class BDI. The superconductivity in 
the artificial heavy-fermion superlattice CeCoIns/YbCoIns and bilayer interface SrTiOs/LaAlOa 
is discussed. 


I. INTRODUCTION 

Parity is an essential quantum number of quantum 
phases unless inversion symmetry is broken. Classifi¬ 
cation of superconducting states is based on the parity 
of the order parameter i According to the conventional 
understanding^ even-parity superconductivity is realized 
by the condensation of spin-singlet Cooper pairs, while 
odd-parity superconductivity is induced by spin-triplet 
Cooper pairs, because of the anticommutation relation 
of fermions. Even-parity superconductivity has been ob¬ 
served in a variety of materials, e.g., archetypal strongly 
correlated electron systems such as high-T c cuprate su¬ 
perconductors (SCs)^£ as well as conventional SCs sta¬ 
bilized by electron-phonon coupling. On the other hand, 
only a few materials are considered as possible hosts of 
odd-parity superconductivity. This is probably because 
the conditions for spin-triplet pairing are unfavorable in 
most materials. Since electron-phonon coupling mostly 
stabilizes spin-singlet s-wave superconductivity, strong 
electron correlation is required for the glues of spin-triplet 
Cooper pairs. However, d -wave superconductivity is sta¬ 
ble in most strongly correlated electron systems 4 

Odd-parity superconductivity has been attracting at¬ 
tention because of its multicomponent order parame¬ 
ters that give rise to multiple superconducting/superfluid 
phases and intriguing phenomena related to sponta¬ 
neous symmetry breaking*^ - — Furthermore, a great 
deal of attention has recently been paid to odd-parity 
SCs because they are candidates for topological su¬ 
perconductivity*!^— However, only Sr 2 Ru 04 ^^ and 
some uranium-based heavy-fermion compounds such as 
UPta&i, UGe^, URhGe, and UCoGei£d£ show strong 
evidence for spin-triplet pairing. 

Recent theoretical studies have presented another way 
to stabilize the odd-parity superconducting state. It has 
been shown that odd-parity superconductivity may oc¬ 
cur through spin-singlet Cooper pairs in crystals lacking 
local inversion symmetry. Such locally noncentrosym- 


metric crystals have a sublattice degree of freedom in 
electronic structures, allowing the odd-parity spin-singlet 
superconducting state . 19 ’ 20 Although this state is not al¬ 
lowed in the absence of spin-orbit coupling according to 
the BCS theory, such an exotic superconducting state 
may be stabilized by the sublattice-dependent spin-orbit 
coupling arising from the relativistic effect.— It has been 
shown that a long-range Coulomb interaction stabilizes 
odd-parity superconductivity in combination with spin- 
orbit coupling —’ 22 On the other hand, two of the au¬ 
thors have shown that the odd-parity spin-singlet super¬ 
conducting state is stabilized by spin-orbit coupling and 
the paramagnetic effect without relying on the particu¬ 
lar electron correlation effect— Therefore, conventional 
electron-phonon coupling or antiferromagnetic spin fluc¬ 
tuation leading to spin-singlet Cooper pairing may induce 
odd-parity superconductivity when both spin-orbit cou¬ 
pling and the paramagnetic effect play important roles. 

In a previous study, we focused on two-dimensional 
(2D) multilayer SCs, in which global inversion symmetry 
is preserved but some of the layers lack local inversion 
symmetry. Then, by applying a magnetic field along the 
c axis, the order parameter of spin-singlet superconduc¬ 
tivity changes sign across the center layer, as shown in 
Fig. 1(a)— The order parameter spatially modulated in 
the atomic length scale ensures the odd parity of super¬ 
conductivity. Such a superconducting state is called the 
pair-density-wave (PDW) state. Interestingly, the PDW 
state is classified into topological crystalline supercon¬ 
ductivity protected by mirror reflection symmetry when 
the number of superconducting layers is odd>2 4 A non¬ 
trivial topological invariant in the symmetry class D, the 
mirror Chern number, ensures the appearance of the Ma- 
jorana edge mode at the edge. A promising candidate for 
realizing such a topological superconducting state is the 
recently grown artificial superlattices CeCoIn 5 /YbCoIn 5 
composed of the quasi-2D heavy-fermion SC CeCoIns 
and the conventional metal YbCoIns— The supercon¬ 
ductivity occurs in CeCoIn ,5 multilayers, and YbCoIns 
plays the role of spacer layers. Thus, the superlattice 


2 


(b) Tricolor superlattice 

oooo'oooo 

oooooooo 

oooooooo 

oooooooo 

•••••••• 

oooooooo 

•••••••• 

oooooooo 

•••••••• 

oooooooo 

oooooooo 

oooooooo 

oooooooo 


FIG. 1. (Color online) Illustration of superlattice structures 
studied in this paper. We show examples of the trilayer sys¬ 
tem (M = 3). (a) Binary superlattice that has been realized 
in CeCoIns/YbCoIns. Closed (blue) and open (red) circles 
show the superconducting layers and spacer layers, respec¬ 
tively. The layer-dependent order parameters in the BCS and 
PDW states are described on the left, (c) Tricolor superlat¬ 
tice. Thick (red) and thin (green) open circles show the two 
kinds of spacer layers. 


is regarded as a 2D multilayer SC when the number 
of YbCoIn 5 layers is large. The strong spin-orbit cou¬ 
pling and the large paramagnetic effecl— ^ in CeCoIn 5 
are favorable for the topological odd-parity supercon¬ 
ductivity discussed above. Indeed, the effects of layer- 
dependent Rashba spin-orbit coupling (RSOC) on the 
superconducting state of CeCoIn 5 /YbCoIn 5 have been 
observed i 28 ! 29 

On the other hand, assuming a large Maki parameter 
«m = V%H° 2 h /H^ 2 , Ref. 23 neglected the orbital effect 
that competes with the paramagnetic effect. H °£ b and 
Hf 2 are fictitious upper critical fields determined by the 
orbital effect (orbital limit) and by the paramagnetic ef¬ 
fect (paramagnetic limit), respectively. Since the Maki 
parameter of bulk CeCoIn 5 is moderate («m ~ 3) for 
the magnetic field along the c axis^S the orbital effect 
may suppress the PDW state. A simple way to reduce 
the orbital effect is to apply a magnetic field along the 
2D conducting plane. However, when the orbital effect is 
completely neglected, the in-plane magnetic field induces 
the complex stripe (CS) stated and the PDW state is 
not stabilized. In this paper, we show that the odd-parity 
PDW state is stabilized by switching on a weak orbital 
effect. 

Furthermore we find that the PDW state is a topolog¬ 
ical crystalline SC protected by magnetic mirror symme¬ 
try. Because mirror reflection symmetry with respect to 
the ab plane is broken by the in-plane magnetic field, the 
mirror Chern number is no longer a topological invariant. 
However, magnetic mirror symmetry is preserved when 
we apply the magnetic field along the a axis or the b axis. 
Using this symmetry, we define the topological invariant, 
i.e., the one-dimensional winding number in the symme¬ 
try class BDI. We will show a non-trivial winding number 


and the resulting Majorana edge state. 

The organization of this paper is as follows. In Sec. II 
we introduce the model for multilayer SCs possessing the 
layer-dependent RSOC. Various superconducting states 
obtained by solving the linearized mean-field equation 
are illustrated. We study the binary superlattices in 
Sec. Ill and discuss the superconducting state in the ar¬ 
tificial superlattice CeCoIns/YbCoIns. The tricolor su¬ 
perlattice in Fig. 1(b) is also investigated in Sec. IV. In 
Secs. Ill and IV, it is shown that the odd-parity PDW 
state is stabilized through the RSOC and the orbital ef¬ 
fect. In Sec. V, the topologically non-trivial properties 
of the PDW state are clarified, and the Majorana edge 
state is demonstrated. A brief summary and discussions 
are provided in Sec. VI. 

II. MODEL AND FORMULATION 
A. Model 

First, we introduce the model for 2D multilayer SCs. 
By simply neglecting the spacer layers, the binary and 
tricolor superlattices in Fig. 1 are described by the multi¬ 
layer model. By taking into account the layer-dependent 
RSOC, the orbital effect, and Zeeman coupling, the 
Hamiltonian is described as 

E ^ k +Pm)c{ srn c k sm + t± E C ksm C ksm' 

k,s,m fc, s,(m,m') 

+ V(fc, k )Cfc + ^ m cL fc _| m c_fc'_4. m Cfc^ m 

k.k' ,q,m 

— Mb H ■ cr ss ic\. srn Cks l m 

k,s,s' ,m 

H” ^ ^ Pm) * &ss' ^ksm^ks'mi (1) 

k,s,s' ,m 

where k, s, and m (= 1,..., M) are indexes of momentum, 
spin, and layer, respectively. The number of supercon¬ 
ducting layers is M. 

The first term is the energy dispersion in the single¬ 
layer limit. We adopt the nearest neighbor hopping term 
in the square lattice for simplicity, 

£(fc) = — 2t(cos k x a + cos k y a) — p, ( 2 ) 

where a is the lattice constant. The orbital effect in¬ 
duced by the applied magnetic field is taken into ac¬ 
count through the Peierls phase. When we consider the 
magnetic field along the [100]-axis, H = Hx, we can 
choose the vector potential A = (0, — Hz,0). Then, 
the orbital effect leads to the layer-dependent shift of 
momentum, k —> k + j-A. Thus, we obtain p m = 
j^Hd [m — ( M + l)/2] y , where d is the lattice spacing be¬ 
tween the nearest neighbor superconducting layers. For 
binary superlattices, d = c with c being the lattice con¬ 
stant along the c-axis. Later we adopt the lattice con¬ 
stant of CeCoIns, a = 4.6 A and c = 7.5 A since we 
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focus on the artificial superlattice CeCoIns/YbCoIns^ 
The second term of Eq. m is the inter-layer hopping 
term. Since we consider the heterostructures composed 
of quasi-2D SCs, the inter-layer hopping t± is assumed 
to be much smaller than the in-plane hopping, t± <C t. 

The third term describes the pairing interaction. We 
assume s-wave superconductivity for simplicity, and thus, 


V(k 1 k’) = -V a . 


( 3 ) 


We have confirmed that qualitatively the same results are 
obtained for d- wave superconductivity. As we will show 
layer, a spatially non-uniform superconducting state may 
be stabilized. Thus, we take into account the finite 
center-of-mass momentum of Cooper pairs q , and k± = 
k ± q/2. The fourth term is the Zeeman coupling term 
giving rise to the paramagnetic effect on the supercon¬ 
ducting state. 

We show that exotic superconducting states are sta¬ 
bilized by the layer-dependent RSOC represented in the 
last term of Eq. (1), which arises from the local viola¬ 
tion of inversion symmetry^ Ensured by the global in¬ 
version symmetry in the crystal structure, the RSOC is 
odd with respect to the mirror reflection on the center 
layer, and thus ctM+i-m = —c* m - For example, the layer- 
dependent coupling constant is ( 011 , 0 : 2 ) = (a,—a) for 
bilayers and ( 01 , 02 , 0 : 3 ) = (o,0, — a) for trilayers. We 
assume a g-vector characterizing the RSOC^ g(k) = 
(— sin k y a, sinfc x a, 0 ), so as to satisfy the periodicity in 
momentum space. 

For bilayers, a similar model has been investigated 
by noticing the twin boundary of noncentrosymmetric 
SCsi22r— Then, the intriguing superconducting phase 
with broken time-reversal symmetry was investigated by 
assuming comparable pairing interactions for the spin- 
singlet Cooper pairs and spin-triplet onesi^^r— We avoid 
such fine-tuning of pairing interactions here and consider 
the dominantly spin-singlet pairing state, which is real¬ 
ized in most SCs. Aoyama et. al. studied the magneto¬ 
electric effect on the upper and lower critical magnetic 
fields^ but their Ginzburg-Landau model does not ap¬ 
propriately take into account the paramagnetic effect in 
the high magnetic field region, and therefore, the PDW 
state focused in this paper is not obtained. 

We assume small inter-layer hopping t±/t = 0.1 and 
a moderate RSOC a/t = 0.3, unless explicitly men¬ 
tioned otherwise. As shown by previous studies ; 20 ’ 23 ! 31 
exotic superconducting states may be stabilized when 
M/t± > 1. Thus, we assume |a|/tj_ > 1 throughout this 
paper. This condition may be satisfied in heterostruc¬ 
tures of quasi-2D compounds. The chemical potential 
is g/t = 2 in Sec III and IV while it is g/t = —2 in 
Sec V. We choose the pairing interaction V s /t = 1.3 or 
V s /t = 1.5. We confirmed that the following results are 
almost independent of the choice of V s /t. 


B. Linearized mean-field theory 


We study the superconducting state by means of the 
linearized mean-field theory. Although we have to fully 
solve the Bogoliubov-de Gennes (BdG) equation in or¬ 
der to obtain the superconducting phase diagram, we can 
clarify the superconducting state near the transition tem¬ 
perature by linearizing the BdG equation while avoiding 
the numerical limitations of the full BdG equation. The 
linearized BdG equation is formulated by calculating the 
superconducting susceptibility, 

X s *m'(<l)= f dTe^{B qm {T)Bl m ,( 0)), (4) 

Jo 

where q = (q,ifl n ) and Q n = 2nnkBT is the boson Mat- 
subara frequency. The annihilation operator of Cooper 
pairs is introduced as 


B 


qm — 


'y ' Ck+q\mC—k\,m-> 

k 


( 5 ) 


and B qm (r) = e nT B qrn e HT . 

The superconducting susceptibility is obtained by us¬ 
ing the T-matrix approximation, 


X sc (q) = 


X°(g) 

i-14x°(g)’ 


( 6 ) 


where y sc = (Xmm') the Af x M matrix. The irre¬ 
ducible susceptibility is calculated by 

Xmm> (<?) = ~ ^l) 

k,l 

-Gw ( fc + 9’ *<*>«) G w (“ fc ! in n - m)\ , (7) 

where G^ m , (k,iu>i) is the non-interacting Green func¬ 
tion and uji = (21 + l^k^T is the fermion Matsubara 
frequency. 

Superconducting transition occurs at the temperature 
where X sc (q) diverges. Thus, the criterion of the super¬ 
conducting instability is obtained as the largest eigen¬ 
value of I4x°(<z) is unity. We obtain the layer-dependent 
order parameter A m (r) = A,„e lq r from the eigenvector, 
(Ai, A 2 , ...Am) t . Because the global inversion symme¬ 
try is preserved in our model, the eigenvalues are equiv¬ 
alent between the momentum ±q. Hence, the single-g 
state or the double-g state may be stabilized when the 
center-of-mass momentum of Cooper pairs is finite. It is 
expected that the double-q state is stable in our model 
because the order parameter almost disappears in one of 
the outermost layers in the single-g state with a small 
condensation energy. For instance, Am "C Ai for mo¬ 
mentum q while Ai -C Am for the opposite momentum 
— q. In the double-g state, the order parameter is de¬ 
scribed as A m (r) = A^e lq ' r + AmV' q ' r , where A^ 
is the eigenvector of V s x°(q) for the momentum ±q, re¬ 
spectively. We confirmed that the bosonic Matsubara 
frequency is always zero, il n = 0. 
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C. Superconducting states 

In this subsection, we classify the solution of the lin¬ 
earized BdG equation. As we will show later, various su¬ 
perconducting states are stabilized in our model. They 
are illustrated for the bilayer system in Fig. [2] In this sub¬ 
section we discuss the bilayer system for simplicity, since 
the extension to more-than-two-layer systems is straight¬ 
forward. 


(a) BCS state 

(b) Vortex state 

(c) CS state 

(d) PDW state 
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FIG. 2. (Color online) Illustration of superconducting states 
in the bilayer system, (a) BCS state, (b) vortex state, (c) CS 
state, and (d) PDW state. Thick bars show the superconduct¬ 
ing layers. Layer-dependent order parameters are described 
in the figures. The spatial dependence of the superconduct¬ 
ing gap is illustrated by thin lines on top of the figures. The 
vortex in (b) and anti-vortex in (c) are shown by arrows. 


The uniform superconducting state [A m (r) = A] is 
stable at zero magnetic field as expected from conven¬ 
tional BCS theory. Thus, we call the uniform state the 
“BCS state” [Fig. [U[a)]. On the other hand, a variety 
of spatially non-uniform states may be stabilized in the 
magnetic field. First, the orbital effect induces the vortex 
state illustrated in Fig. 0(b). When the quantum vortices 
penetrate inside multilayers, the order parameter is de¬ 
scribed as 

Ax(r) = A 0 (e _iq ' r + <5e iqr ), (8) 

A 2 (r) = Ao(5e- jqr + e iqr ), (9) 

where |<S| < 1, and q = q(z x H) with H = H/\H\ 
and q > 0. Second, the layer-dependent RSOC stabilizes 
the CS state through the paramagnetic effect^ The CS 
state is also described by Eqs. © and m, however, the 
sign of the center-of-mass momentum q depends on the 
band structure and the sign of RSOC. For our choice of 
parameters, the RSOC favors q < 0 when a > 0 while 
q > 0 when a < 0. Thus, the CS state is regarded 
as an anti-vortex state [Fig. O c )] and is distinguished 
from the vortex state, when a > 0. Then, the spin- 
orbit coupling competes with the orbital effect, and gives 
rise to an intriguing superconducting phase diagram, as 
we show later. We focus on this case in the following, 
although the other case, a < 0, is briefly discussed. 

Finally, Fig. [2](d) illustrates the PDW state, where 
[Ai(r), A 2 (r)] = (A, — A)32 The order parameter is uni¬ 
form in the 2D conducting plane, but it changes sign be¬ 
tween layers. As we mentioned before, the PDW state is 
an odd-parity superconducting state, although the super¬ 
conductivity is induced by the spin-singlet s-wave Cooper 


pairs. Although we have shown that the PDW state 
is stabilized in the c-axis magnetic field near the Pauli 
limit^ in this paper we show that the PDW state is also 
stabilized in the in-plane magnetic field when the RSOC 
competes with the orbital effect. 



Bilayer 

Trilayer 

BCS 

Ai(r) = A 

A 2 (r) = A 

Ar(r) = A 

A 2 (r) = A' 

A 3 (r) = A 

PDW 

Ai(r) = A 

A 2 (r) = -A 

Ai(r) = A 

A 2 (r) = 0 

A 3 (r) = —A 

FFLO 

Ai(r) = Acos(q ■ r) 

A 2 (r) = Acos(q ■ r) 

Ai(r) = Acos(q ■ r) 

A 2 (r) = A' cos(q • r) 

A 3 (r) = Acos(q ■ r) 

Vortex 

At(r) = A(e _iqr + 5e iqr ) 
A 2 (r) = A((5e _iqr + e iqr ) 

At(r) = A(e _iqr + <5e iq ' r ) 
A 2 (r) = A' cos(q • r) 

A 3 (r) = A(5e iqr +e iqr ) 

CS 

Ai(r) = A(e iqr + <5e _iqr ) 
A 2 (r) = A(5e iqr + e _iqr ) 

Ai(r) = A(e iqr + 5e _iq ' r ) 
A 2 (r) = A' cos(q • r) 

A 3 (r) = A(5e iqr + e _iqr ) 


TABLE I. Layer-dependent order parameter of the BCS, 
PDW, FFLO, vortex, and CS states in bilayers and trilay¬ 
ers. 


The BCS, PDW, vortex, and CS states are distin¬ 
guished from each other in more-than-two-layer systems 
too. In Table I, we summarize the order parameter 
of these states in trilayers as well as in bilayers. The 
order parameter of the so-called Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) stat e 37 ’ 38 is also shown for a com¬ 
parison. 


III. BINARY SUPERLATTICE 
A. Bilayer system 

In this section, we study the binary superlat¬ 
tices (see Fig. [2(a)) which have been fabricated in 
CeCoIn 5 /YbCoIn 5 .^ 

We begin with the simplest case, namely, the bilayer 
system (M = 2) illustrated in Fig. [2] Then, the strength 
of the orbital effect is controlled by the Fermi energy 
Ep , which is proportional to the in-plane hopping t. 
The orbital effect is estimated from the dimensionless 
quantity U£c/ < f’o with <f>o = being the flux quan¬ 
tum. Since the coherence length of superconductivity 
is £ ~ hvY/k^,T c ~ [E-p/k^T c )a, the orbital effect is en¬ 
hanced by increasing the Fermi energy. On the other 
hand, the paramagnetic effect of the magnetic field is es¬ 
timated from another dimensionless quantity /luH / kuT c . 
that is independent of the Fermi energy. Thus, the orbital 
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effect (paramagnetic effect and RSOC) plays an impor¬ 
tant role in the superconducting state for large (small) 
in-plane hopping t. Note that the RSOC induces the CS 
state through the paramagnetic effect. 


(a) 
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FIG. 3. (Color online) (a) Transition temperatures of vari¬ 
ous superconducting states in the bilayer system. We assume 
a small in-plane hopping, t = 20 meV. Thin (red), moderate 
(blue), and thick (green) lines show the T c of BCS, CS, and 
PDW states, respectively. The highest transition temperature 
is observable and drawn by the solid line, while the dashed 
and dot-dashed lines show the fictitious transition temper¬ 
atures. (b) Magnetic field dependence of the center-of-mass 
momentum |q|. (c) Layer dependence of the order parameter. 
We show A 2 /A 1 for the BCS state and S for the CS state. In 
this subsection, we choose the pairing interaction V s /t = 1.5, 
which gives the transition temperature AnTco = 0.01241 in the 
absence of the spin-orbit coupling and magnetic field. The 
temperature and magnetic field are scaled by T c o in the fig¬ 
ures. 


When we assume small in-plane hopping t = 20 meV 
consistent with the heavy effective mass of CeCoIns, the 
orbital effect is negligible. Indeed, we obtain the phase 
diagram in Fig. 0a) which resembles to the result in the 
paramagnetic limit 41 The CS state is stable in the high 
magnetic field region, but the PDW state is not stabi¬ 
lized. As increasing the magnetic field, the center-of-mass 
momentum of Cooper pairs gradually increases through 
the BCS-CS phase transition (Fig. [3](b)) , indicating the 
second order phase transition. Figure 0c) shows that 
Ai = A 2 in the BCS state while <5 in Eqs. (0) and (0) 
decreases with increasing the magnetic field in the CS 
state. These behaviors are consistent with the previous 
study on the same model^i where the orbital effect is sim¬ 
ply neglected. Thus, the previous study that was focused 
on the heavy-fermion superlattice CeCoIns/YbCoIns is 


justified. 


(a) 



FIG. 4. (Color online) Transition temperatures (a), center-of- 
mass momentum (b), and layer-dependent order parameter 
(c) for a large in-plane hopping, t = 200 meV. The other 
parameters are the same as those in Fig. [3] The moderate 
(purple) line in (a) shows the transition temperature of the 
vortex state. 


On the other hand, the orbital effect significantly af¬ 
fects the superconducting state for large in-plane hop¬ 
ping, t = 200 meV. Figure 0a) shows that the vortex 
state is stable in the high magnetic field region, as ex¬ 
pected. The BCS-vortex phase transition is second or¬ 
der as indicated by the continuous change of the center- 
of-mass momentum [Fig. 0b)] and S [Fig. 0c)]. Al¬ 
though the RSOC and the paramagnetic effect play less 
important roles than the orbital effect, they induce a 
characteristic magnetic field dependence in the center- 
of-mass momentum around the Pauli-Chandrasekhar- 
Clogston limit ijl-qH/ k^,T c0 = 1.25. The RSOC and para¬ 
magnetic effect suppress the orbital effect, and thus de¬ 
crease the center-of-mass momentum above the Pauli- 
Chandrasekhar-Clogston limit. 

A main result of our study is obtained when the RSOC 
competes with the orbital effect. Such a situation is re¬ 
alized for moderate in-plane hopping t = 80 meV. Then, 
the PDW state is stabilized in the high magnetic held 
region, as shown in Fig. 0 Note that the PDW state is 
induced neither by the orbital effect nor by the paramag¬ 
netic effect and RSOC. The PDW state is stable due to 
a balance of these effects. We explain this mechanism in 
details here. The order parameter of the CS and vortex 
states is described by Eqs. © and ©. and these two 
states are differentiated by the sign of q. The positive q 
(vortex state) is favored by the orbital effect, although 
the negative q (CS state) is induced by the RSOC. Thus, 
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PDW state is stabilized. 

We take into account the RSOC on the outermost lay¬ 
ers, (au, q;m) = (a,—a), while the RSOC on the other 
layers is neglected for simplicity. This is a reasonable as¬ 
sumption for the layer dependence of RSOC, because the 
spin-orbit coupling is determined by the local environ¬ 
ment of atoms^ and thus the outermost layers contain 
the largest spin-orbit coupling. 



FIG. 6. (Color online) Transition temperatures of the BCS, 
CS, and PDW states in the trilayer system (M = 3) for small 
in-plane hopping, t = 21 meV. In this subsection, we assume 
V s /t = 1.3, which gives rise to the transition temperature 
at zero magnetic field, T c o = 0.00487f/fcB, in the absence of 
RSOC. 


FIG. 5. (Color online) Transition temperatures (a) and layer 
dependence of the order parameter (b) for a moderate in-plane 
hopping, t = 80 meV. The other parameters are the same as 
those in Fig. [3] It is shown that the PDW state is stable in 
the high magnetic field region. 


q ~ 0 when these two effects are in balance. Then, the 
inter-layer Josephson coupling stabilizes the uniform su¬ 
perconducting state along the conducting plane, and thus 
q = 0. The 7r phase difference between layers is favored 
so that the paramagnetic depairing effect is avoided as in 
the c-axis magnetic fields In this way, the PDW state is 
stabilized by the orbital effect, the paramagnetic effect, 
RSOC, and inter-layer coupling. 

B. More-than-two-layer system 

Although it is hard to experimentally control the Fermi 
energy, the number of superconducting layers M can be 
tuned by using the artificial superlattice j 25 ! 28 ’ 29 Thus, 
we may be able to control the orbital effect by tuning 
M. Since the shift of momentum on the outermost lay¬ 
ers, |pi| = \pm\ = eHc(M — l)/2 h, increases with M, 
the orbital effect is enhanced by increasing the number 
of superconducting layers. This is reasonable because 
vortices easily penetrate inside of thick SCs. We demon¬ 
strate here that the competing region of the RSOC and 
the orbital effect is realized by tuning M, and then the 


We fix the in-plane hopping, t = 21 meV, and the 
pairing interaction, V s /t = 1.3. Then, the orbital effect is 
negligible in the bilayer system as in Fig. [3] The trilayer 
system shows the phase diagram (Fig. [G]) similar to Fig.(3j 
and thus the trilayer system is still close to the Pauli 
limit. 



FIG. 7. (Color online) Transition temperatures of the BCS, 
CS, and PDW states in the five-layer system (M = 5) for 
small in-plane hopping, t = 21 meV. 


On the other hand, the orbital effect plays an impor¬ 
tant role in the five-layer system (M = 5). Then, the 
orbital effect competes with the RSOC, and therefore, 
the PDW state is stabilized as expected from the results 
in Sec. IIIA. Indeed, Fig. Q] shows that the BCS state, 
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CS state, and PDW state are stabilized in the low, in¬ 
termediate, and high magnetic field regions, respectively. 
Because the CS-PDW transition is first order phase tran¬ 
sition, the upper critical field shows a kink, although 
the kink may be weak in some cases (see Figs. 5(a) and 
10(a)). Generally speaking, a distinct kink appears when 
the transition temperature of the PDW state is small. 
Then, the H c 2 of the CS state is suppressed, while that 
of the PDW state shows an upward curvature. The ob¬ 
servation of the kink will be an experimental test for the 
presence of the PDW state. 



T/Tco 


(a) (b) 




T/T co T/T co 

FIG. 9. (Color online) Transition temperatures of the super¬ 
conducting states in the five-layer system, (a) Pauli limit by 
setting pm = 0. (b) Orbital limit by eliminating the Zeeman 
term in the Hamiltonian [Eq. ©]• The other parameters are 
the same as those in Fig. [7] 


field through the crossover. We see that the paramag¬ 
netic effect appears even in the seven-layer system, where 
the orbital effect is larger than the effect of RSOC. 


C. CeCoIns/YbCoIng 


FIG. 8. (Color online) Transition temperatures of the BCS, 
vortex, and PDW states in the seven-layer system (M = 7) 
for small in-plane hopping, t = 21 meV. 


When we furthermore increase the number of super¬ 
conducting layers, the superconducting state is domi¬ 
nated by the orbital effect, and thus the vortex state 
is stabilized in the high magnetic field region. For exam¬ 
ple, we show the phase diagram of the seven-layer system 
(M = 7) in Fig. [5] It is shown that the PDW state is not 
stabilized. 

We focus here on the five-layer system, and we empha¬ 
size the cooperative role of the RSOC, the orbital effect, 
and the paramagnetic effect on the thermodynamical sta¬ 
bility of the PDW state. FigureEKa) shows the phase di¬ 
agram in the absence of the orbital effect. The CS state 
is more stable than the PDW state, as in the bilayer and 
trilayer systems. On the other hand, the vortex state is 
stable in the high magnetic field region when we neglect 
the paramagnetic effect, as shown in Fig.[9])b). Note that 
the RSOC does not play important roles in the absence 
of the paramagnetic effect. 

It should be noticed that the orbital limit of the upper 
critical field [Fig. [5][b)] is much larger than the param¬ 
agnetic limit [see Fig. [Ufa)], indicating the large Maki 
parameter. This means that the upper critical field of 
five-layer systems is mainly determined from the param¬ 
agnetic effect. In this sense, the PDW state occurs near 
the paramagnetic limit, although a weak orbital effect is 
needed. Stars in Figs. [Gl El [HI and[9](a) show the crossover 
induced by the paramagnetic effect.— Because the cen¬ 
ter layer is not protected against the paramagnetic effect 
by the RSOC, the order parameter in the center layer 
A(m+i )/2 suddenly decreases by increasing the magnetic 


In the previous subsections, we designed the odd-parity 
PDW state using the artificial heterostructures. It has 
been shown that the orbital effect is controlled by the 
number of superconducting layers. Indeed, various su¬ 
perlattices CeCoIns/YbCoIns with M > 2 are supercon¬ 
ducting, and we can tune the number of layers M 1 25 i 28 i 29 
Thus, the artificial superlattice CeCoIns/YbCoIns may 
be a new platform for the odd-parity superconductivity. 
A reasonable parameter t ~ 20 meV leads to the PDW 
state in the five-layer system. 

We comment here on the recent experimental obser¬ 
vations of the paramagnetic effect in the superlattice 
CeCoIns/YbCoIns— Goh et. al. observed the strong 
paramagnetic effect by measuring the field-angle depen¬ 
dence of the upper critical field. They also showed that 
the paramagnetic effect is suppressed in a few-layer sys¬ 
tem M < 3. Their experimental results are consistent 
with our model; the paramagnetic effect is suppressed 
with decreasing M because the superconductivity in sur¬ 
face layers, m = 1 and M, is substantially protected 
against the paramagnetic effect due to the RSOC .—' — 
Indeed, the upper critical field of the trilayer system is 
larger than that of the five-layer system near T = T c 0 
(see Figs. 6 and 7). Note that the upper critical field 
is dominantly determined from the paramagnetic effect 
even when the orbital effect competes with the RSOC 
(see the discussion in Sec. IIIB). 

Considering the consistency between our calculation 
and experiments for CeCoIns/YbCoIns at low magnetic 
fields, it is expected that the PDW state may be realized 
in the artificial superlattice with M ~ 5 at high mag¬ 
netic fields. However, any indication of the presence of 
high-field superconducting phase has not been reported. 
For instance, the kink and the upturn of the upper crit- 













ical field shown in our calculations have not been ob- 
served i 25 ! 28 ! 29 The high-field phase may have been missed 
because the measurement has not been carried out in the 
low-temperature region. On the other hand, the discrep¬ 
ancy may be attributed to the two ingredients that are 
not taken into account in our model. First, it is ex¬ 
pected that disorders suppress the high-field supercon¬ 
ducting phase, as the FFLO state is suppressed^ The 
artificial superlattice CeCoIns/YbCoIns indeed contains 
substantial disorders. Second, the number of spacer lay¬ 
ers TV = 4 ~ 6 may not be large enough to eliminate 
the coupling between superconducting multilayers. Inter- 
multilayer coupling is harmful for the PDW state, and 
thus it should be decreased by increasing the number of 
spacer layers. 


IV. TRICOLOR SUPERLATTICE 

Next, we discuss the tricolor superlattice illustrated 
in Fig. mb). As we have shown in Sec. IIIB, the or¬ 
bital effect is controlled by the spacing between the out¬ 
ermost layers. Thus, the PDW state may be stabilized 
in a tricolor superlattice by intercalating the spacer lay¬ 
ers into the superconducting layers. Then, the spac¬ 
ing of neighboring superconducting layers is multiplied 
to d = (mu + l)c with md being the number of inter¬ 
calated spacer layers (green open circles in Fig. 0b)). 
We assume here that the inter-layer spacing between the 
spacer layer and the superconducting layer is c. A sim¬ 
ilar situation in the bilayer d-doped SrTiO^i has been 
realized, and superconductivity in the bilayer interface 
LaAlOs/SrTiOs/LaAlOa has been studied^. 

Multiplying the inter-layer spacing by (m<j+l) is equiv¬ 
alent to increasing the in-plane hopping to (md+l)t while 
keeping the ratio, t±/t, n/t, a/t , and V s /t. For instance, 
we obtain the same results for the binary superlattice 
with t = 80 meV and for the tricolor superlattice with 
t, = 20 meV and m d = 3. Thus, the PDW state may 
be stabilized in the tricolor bilayer superlattice. How¬ 
ever, in reality, the inter-layer hopping t± between the 
nearest-neighbor superconducting layers is significantly 
decreased by intercalating a spacer layer. For instance, 
we obtain t± ~ t^ s /(E s — E n ) in the presence of a single 
spacer layer (md = 1), where t ns is the hopping integral 
between the superconducting layer and the spacer layer 
and E s — E n is the potential difference between these 
layers. Therefore, the stability of the PDW state in the 
tricolor superlattice should be examined by investigating 
the superconducting state for small t±. 

We study here the tricolor bilayer superlattice for sim¬ 
plicity. The in-plane hopping is set to t = 20 meV by con¬ 
sidering the heavy-fermion superlattice. As we discussed 
above, we obtain the same phase diagram as Fig. [5] when 
Wd = 3 and t±/t = 0.1. On the other hand, Figs. fTOl al 
and 10(b) show the phase diagram for t±/t = 0.05 and 
t±/t = 0.01, respectively. Since the inter-layer Josephson 
coupling decreases with t±/t, the uniform states, namely, 


(a) 



(b) 



FIG. 10. (Color online) Transition temperatures of the BCS, 
CS, vortex, and PDW states in the tricolor superlattice with 
md = 3, t as 20 meV, fj,/t = 2, a/t = 0.3, and V s /t = 1.5. 
(a) t±/t = 0.05 and (b) t±/t = 0.01. The temperature is 
normalized by T c o = 0.0124t/fcB- 


the BCS and PDW states, are suppressed. It is shown 
that the PDW state is stable at high magnetic fields, but 
the transition temperature of the PDW state is signifi¬ 
cantly decreased for t±/t = 0.01. The PDW state will 
be furthermore suppressed by further decreasing t±/t. 
Thus, it may be hard to realize the PDW state in a tri¬ 
color superlattice by intercalating many spacer layers. 

It should be noticed that Fig. Mb) is similar to 
the phase diagram of the binary five-layer superlattice 
(Fig-0- We now understand that the reduced transition 
temperature of the PDW state in the five-layer system is 
due to the reduced Josephson coupling between the out¬ 
ermost layers. In other words, the superconducting inner 
layers play a role of the spacer layers. Indeed, the PDW 
state is mainly induced by the outermost layers where 
the superconductivity is protected against the paramag¬ 
netic effect by the RSOC. For example, we obtain the 
layer-dependent order parameter (Ai, A 2 , A 3 , A 4 , A 5 ) ~ 
A(l, 0.18, 0, —0.18, —1) in the five-layer PDW state. 

Finally, we discuss the superconductivity in the 
bilayer (5-doped SrTiO^ and the bilayer interface 
LaA 103 /SrTiC> 3 /LaA 103 .— The Fermi velocity esti¬ 
mated for a carrier density n ~ 10 14 cm -2 on the ba¬ 
sis of the three-orbital tight binding model - 42 is approx¬ 
imately twice as large as that in our model. Since the 
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lattice constant along the c-axis c = 3.9 A is nearly half 
of CeCoIn 5 , the orbital effect in the SrTi0 3 heterostruc¬ 
tures is comparable to our model. Thus, the bilayer sys¬ 
tem sandwiching three non-superconducting layers may 
be a platform of the odd-parity PDW state. Then, the 
inter-layer coupling between superconducting layers may 
be small, and therefore, the PDW state may appear in 
the low temperature region, as shown in Fig. 10(b). The 
kink in the upper critical field would be a signature of 
the PDW state. 


V. TOPOLOGICAL SUPERCONDUCTIVITY 


Topologically non-trivial insulators and SCs have 
evolved into one of the major research topics of mod¬ 
ern condensed matter physics recently . 43 ^ 44 In particu¬ 
lar, topological superconductivity attracts a great deal 
of attention since the Majorana state satisfying the 
non-Abelian statistics appears at the edge and disloca- 
tions AIT— In addition to the “strong” topological phases 
classified based on the topological Periodic TablerAA 4 , 
theories on the symmetry-protected topological super¬ 
conducting phases (topological crystalline superconduc¬ 
tivity) have developed recentlyASc— In the c-axis mag¬ 
netic field the odd-parity PDW state is a topological 
crystalline SC protected by mirror symmetry ^ 4 The mir¬ 
ror symmetry along the ab plane protects the topologi¬ 
cal invariant, that is mirror Chern number— of symme¬ 
try class D. We obtain the finite mirror Chern num¬ 
ber i/(±z) = =pl, marking the topologically non-trivial 
properties of superconductivity, when the number of su¬ 
perconducting layers M is oddAA On the other hand, the 
mirror Chern number is no longer a topological invariant, 
when the magnetic field is tilted from the c axis. 

In this section we show that the PDW state may be¬ 
long to another kind of topological crystalline SC when 
the magnetic field is applied along the a or 6 axis. We 
demonstrate the topologically non-trivial properties on 
the basis of the BdG Hamiltonian 


'%BdG = + Pm)c l ksrn Ck srn + t_ L E C ksm Cksm ' 


k,s,m 


k,s,(m,m') 


E MB H ■ (Tss'C^mCks'r 


k,s,s' ,m 


j 


+ ^ ^ H - Pm) ’ &ss' ^k srn pks'r 

k,s,s' ,m 

+ o E ^ss'm(fc)Cfc sm cL fcs , m + h.c 


( 10 ) 


where A m (fc) = (A ss / m (fc)) = [i/j m + d m (k) ■ cr]ia y de¬ 
scribes the layer-dependent order parameter of supercon- 
ductivityA Although the purely s-wave superconductivity 
is considered in Sec. Ill and IV, the p- wave component is 
admixed through the layer-dependent RSOC by the local 
violation of inversion symmetryAi The layer dependence 


of order parameter is obtained as 

(V’eV’ 2,'03,'04,'05) = (VW^injO, -Ipout) , (11) 

(d 1 (k),d 2 (k),d 3 (k),d 4 (k),d 5 (k)) 

— (dout, din, din, An, dout) d{k)j ( 12 ) 
in the five-layer PDW state. The BdG Hamiltonian is 
represented in Nambu space 


%BdG = o E C^BdG(&)c, 


-ffBdG(fc)= 


(H 0 {k) A(fc) 

V At(fe) -Ho(-k) 1 


(13) 

(14) 


where H 3 (k) is the Hamiltonian in the normal state and 
A(fc) is the order parameter. 

Although the mirror symmetry with respect to the ab 
plane M a b is broken in the in-plane magnetic field, the 
magnetic mirror symmetry T 1 = TAi ca {TM.bc) is pre¬ 
served in the a-axis ( 6 -axis) magnetic field. For instance, 
in the magnetic field along the a axis the BdG Hamilto¬ 
nian is invariant under magnetic mirror symmetry 

T'H B dG{k X: k y )T'^ = H B dG{—k x ,ky), (15) 

where A i C a = i&y is the mirror reflection operator and 
T = ia y )C is the time-reversal operator with K. be¬ 
ing the complex conjugate operator. Combining with 
the particle-hole symmetry CH B dG(k)C t = —H B dG{—k), 
where C = t x K, and t x is the Pauli matrix in the particle- 
hole space, we can define the mirror chiral symmetry 
TH B dG(k x ,k y )T^ = -H Bd G(k x ,-ky) with V = ~CV = 
t x . Thus, the BdG Hamiltonian satisfies the chiral sym¬ 
metry 

{T,H BdG (k)} = 0 (16) 

at k y = 0 and k y = n/a. The chiral symmetry ensures 
that the one-dimensional winding number 

J dk x TV [q(k)~ l d kx q{k) - (fc ) -1 d kx q^ (fc)] 

(17) 

is a topological invariant^ — when a finite gap is open 
at k y = 0 and k y = n/a. The 2 M x 2 M matrix q(k) is 
obtained by carrying out the unitary transformation 


(18) 

When we regard the magnetic mirror symmetry 
T 7 as pseudo-time-reversal symmetry; 58 » 59 the one- 
dimensional Hamiltonian ( k x ) = H B dG{k Xl T) and 

~ n / a (k x ) = H B dG{k x ,ir/a) belong to the symmetry 
class BDI because T' 2 = +1AAAA Thus, we can define 
the integer topological numbers of the BDI class 

Vq DI = — ( dk x Tr [q(k x ,0)~ 1 dk x q{k x ,0)\ , (19) 

m Jo 

v ^fa = —■ f dk x Tr [q(k x ,Tr/a)~ 1 d kx q(k x ,Tr/a)\ . (20) 
m Jo 
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Indeed, these winding numbers are equivalent to Eq. m, 
namely, = w 0j7r / a . The pseudo-time-reversal sym¬ 

metry considered here has been used for the definition of 
the integer topological number in one-dimensional semi¬ 
conductor nanowires^ and quasi-one-dimensional d -wave 
superconductors.— The magnetic mirror symmetry is the 
physical origin of this “hidden” time-reversal symmetry. 
The difference of two winding numbers, z^ DI — zA?P r , is 
the strong index of 2D topological crystalline SCs pro¬ 
tected by the magnetic mirror symmetry— 



0 

k a 


y 

FIG. 11. (Color online) Superconducting gap at each k y in the 
five-layer PDW state, defined by E m in(k y ) = Minj^ |i?i(fc)|. 
We assume t = 20 meV, t±/t = 0.1, fi/t = —2, a/t = 0.3, 
and pbH/ t = 0.04 is the magnetic field along the a axis. 
The layer-dependent order parameters are chosen as V'out = 
0.02 and ipi n = 0.0036 consistent with the results in Sec. IIIB 
(Fig. 7). We take into account a small p- wave component d = 
dout = di n = d[n = 0.004 induced by the RSOC— (thick solid 
line), while we obtain the dashed line for d = 0. We also show 
the gap for d = 0.004 and psH/t = 0.001 for a comparison 
(thin solid line). Bulk gap opens in the low magnetic field 
region, but it closes in the high magnetic field region where 
the PDW state is stable. The superconducting gap at k y = 0 
is finite even in the high magnetic field region. 


We now discuss the superconducting gap. Figure |TT] 
shows the gap of the single-particle excitation spectra in 
the five-layer PDW state for each k y , which is defined 
as E min (k y ) = Min ii fcJ.E i (fe)|, with Ei(k) being eigen¬ 
values of the BdG Hamiltonian HedG^)- We here as¬ 
sume p/t = —2 so that the Fermi surface encloses the T 
point (fc = 0). Since there is no Fermi surface along 
k y = 7r/a, the winding number is trivial, zA^ 1 = 0. 

Therefore, we focus on z/® DI . The superconducting gap 
is finite at k y = 0, ensuring the topological protection of 
the winding number z/® DI . At low magnetic fields, the 
superconducting gap is finite in the whole Brillouin zone 
(thin solid line in Fig. 11), and thus, z/® DI is the strong 
topological index. Although the gap at finite k y is closed 
at high magnetic fields owing to the paramagnetic effect, 
the winding number is regarded as a topological number 
of an effective one-dimensional Hamiltonian. 

As we show in Fig. fTST a). z/® DI discretely changes upon 
increasing the p -wave component in the order parameter, 


(a) 



-0.005 0 0.005 


(b) 



FIG. 12. (Color online) (a) Winding number and (b) super¬ 
conducting gap at k y = 0 as a function of the p-wave com¬ 
ponent d = d 0 ut = din = d[n . The other parameters are the 
same as those in Fig. [TT| 


d = d out = c?i n = d[ n . The superconducting gap is closed 
for special values of d where the winding number jumps 
(Fig. [Tiil' b)). We here ignore the layer dependence of 
the p -wave component for simplicity. This assumption 
has been justified by the BdG equation, which shows the 
nearly layer-independent p -wave component in the PDW 
state- 

Switching on a small p -wave component \d\ > O.OS'i/’out, 
we obtain a finite winding number indicating the topo¬ 
logically non-trivial properties [see Fig. [T2T a)]. This is 
in sharp contrast to the 2D Rashba SC, where the p- 
wave component overwhelming the s-wave component is 
required for topological superconductivity—^ This con¬ 
dition is hardly realized in real materials. On the other 
hand, a small p -wave component induced by the RSOC 
causes the PDW state to be topologically non-trivial, be¬ 
cause the s-wave component of the order parameter is 
small on inner layers when M > 3. In Fig. 1131 we show 
the winding number of the bilayer, trilayer, and four- 
layer systems in the PDW state. Although the bilayer 
system (M = 2) is trivial, the topologically non-trivial 
superconducting state is induced by a small p -wave com¬ 
ponent for M > 3. According to the random-phase- 
approximation (RPA) analysis of the three-dimensional 
Rashba-Hubbard model, the induced spin-triplet compo- 
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nent is approximately 20% of the spin-singlet component 
for a moderate RSOCi 39 ’ 61 Thus, the p-wave component 
is likely to be large enough to realize the topological crys¬ 
talline superconductivity protected by magnetic mirror 
symmetry for M > 3. Note that these conditions are dif¬ 
ferent from those for the topological superconductivity in 
the magnetic field along the c axis. The PDW state is a 
topological crystalline superconductor in the c-axis mag¬ 
netic field when the number of superconducting layers M 
is odd* 24 Then, the p -wave component is not required. 


Q 

ffl o 


-0.005 0 0.005 

d 


4 - 4 Bilayer 

o- -o Trilayer 
—I x — x Four-layer 


: i 




FIG. 13. (Color online) Winding number of the bilayer, 
trilayer, and four-layer PDW state. We assume a layer- 
independent p-wave component d = d ou t = dm as in 
Fig. 1121 The s-wave component has a layer-dependence 
{ipi,ip 2 ) = (Vwt, -r/w), = (Vw, 0, -i/w), and 

(-01, -02, "04) = {ipout,ipin,-ipin,-ipout)- The parameters 

are the same as those in Fig. [TT] 


A non-trivial winding number may ensure the pres¬ 
ence of the Majorana edge state according to the bulk- 
edge correspondence. When the magnetic field is applied 
along the a axis, the magnetic mirror symmetry TA4 ca 
is preserved at the edge perpendicular to the [100]-axis 
[(100) edge]. Therefore, the winding number protected 
by this symmetry corresponds to the number of zero- 
energy edge states according to the index theorem^ 
Indeed, we show the Majorana edge states in Fig. [lU 
The trilayer PDW state with a large superconducting 
gap is considered for simplicity of numerical calculation. 
The energy spectrum is calculated in the open bound¬ 
ary condition along the a axis. The layer-dependent 
order parameters are (ipi,ip 2 ,ip 3 ) = (1,0, — l)?/w and 
(di(fc),d 2 (fe),d 3 (fc)) = (1,1,1 )dg(k). We see the two 
Majorana modes around k y = 0 (Fig. fldT a)). when we 
assume a small p- wave component leading to the wind¬ 
ing number r , ® DI = 2. The Majorana states have a 
linear dispersion since the chiral symmetry defined in 
Eq. (fT6l) is not preserved at 0 < \k y \ < it / a. Because an¬ 
other pseudo-time-reversal symmetry, T"H^dG{k)T"' = 
-f^BdG(—fc) with T" = TM a b , is preserved, the two Ma¬ 
jorana states form “Kramers pairs”. 

Figure [Til a) also shows the zero-energy flat band at 
|fcj,a| = 1.35-1.7. This mode is specified by another wind¬ 
ing number u} , k protected by the pseudo-time-reversal 


(a) 



FIG. 14. (Color online) Energy spectra in the trilayer PDW 
state with open boundaries along the a axis. We consider the 
magnetic field along the a-axis H/t = 0.3 (a) and the b 
axis p,BH/t = 0.1 (b). We assume t±/t = 0.1, p/t = —2, 
a/t = 0.3, and choose the s-wave order parameter ip ou t = 0.5 
and the p-wave order parameter d = 0.2. The orbital effect is 
neglected for simplicity. 


symmetry T" ■ The winding number ui' k is defined at all 
k y , and we obtain w' k = — 1 at k y where the flat band 
appears. Hence, the zero-energy flat band does not have 
any degeneracy. 

Finally, we comment on the anisotropic response to the 
external magnetic field. The magnetic mirror symmetry 
TAica is broken when we apply the magnetic field along 
the b axis. Then, the zero-energy Majorana states dis¬ 
appear at the (100) edge, as expected [Fig. |T4Kb)] . This 
field angle dependence is attributed to the Ising character 
of the Majorana state. Similarly, the Majorana mode ap¬ 
pears (disappears) at the (010) edge in the magnetic field 
along the b axis (a axis), because the mirror symmetry 
along the be plane AAbc is preserved at the edge. 


VI. SUMMARY AND DISCUSSION 

In this paper we studied 2D multilayer SCs influenced 
by the layer-dependent RSOC. We showed that the odd- 
parity PDW state is stabilized by the competing spin- 
orbit coupling and the orbital effect in the magnetic held 
along the 2D conducting plane. We also showed that 


































12 


the PDW state is a topological crystalline SC protected 
by the magnetic mirror symmetry when a small p-wave 
component is induced by the RSOC. The Majorana state 
has been demonstrated at the (100) edge [(010) edge] in 
the magnetic field along the a axis (6 axis). 

Our finding paves the way toward realizing odd-parity 
superconductivity without a considerable pairing inter¬ 
action in the spin-triplet channel. Although spin-triplet 
superconductivity is hardly stabilized in most SCs except 
for a few exceptions, our proposal provides an alternative 
way to create odd-parity SC by using the sublattice de¬ 
gree of freedom. 

Indeed, recent developments in the technology of ar¬ 
tificial heterostructures may enable the design of the 
odd-parity PDW state. Superconducting 2D electron 
systems have been fabricated in the oxide interfaces 
SrTiOa/LaAlO^ and SrTiOs/LaTiOa,— gate-tuned 
SrTiOyli and MoS 2 r^ and the heavy-fermion superlat¬ 
tice CeCoIn 5 /YbCoIn 5 i2£ It has been reported that in¬ 
terfacial (intrinsic) spin-orbit coupling significantly af¬ 


fects the superconducting state in SrTiC >3 heterostruc- 
tures^^H and CeCoIn 5 /YbCoIng2&2£ (MoS^^). 
Furthermore, the multilayer structure has been arti¬ 
ficially controlled in CeCoIns/YbCoIn^ and J-doped 
SrTi 034 i Thus, we expect that odd-parity topological 
superconductivity will be created in these systems by 
tuning the multilayer structure and the magnetic field. 
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